Quantum theory of dynamic multiple light scattering 
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We formulate a quantum theory of dynamic multiple light scattering in fluctuating disordered 
media and calculate the fluctuation and the autocorrelation function of photon number operator for 
light transmitted through a disordered slab. The effect of disorder on the information capacity of 
a quantum communication channel operating in a disordered environment is estimated and the use 
of squeezed light in diffusing-wave spectroscopy is discussed. 



Many disordered media (suspensions, emulsions, 
foams, semiconductor powders, biological tissues, etc.) 
scatter light diffusely and appear opaque. Statistical 
properties of light scattered from disordered media are 
of fundamental importance for such applications as wire- 
less communications [1, .2, 3] and different types of diffuse 
optical spectroscopies (diffuse transmission spectroscopy 
diffusing-wave spectroscopy @, etc.). With a few 
notable exceptions [3, 0], most of theoretical studies in 
this field of research have been conducted by treating 
light as a classical wave , without accounting for 

the quantum nature of light. Meanwhile, recent exper- 
iments demonstrate that quantum properties of light in 
disordered media (such as quantum noise 11] and photon 
statistics 12]) can be of interest from both fundamental 
and applied points of view. Moreover, the progress in 
generation of light with nonclassical properties (single 
photons and light in number states [13|, squeezed light 
[3], entangled photons [l^]) suggests that sources of non- 
classical light could soon become accessible to a wide 
scientific community. This calls for a better theoretical 
understanding of quantum aspects of multiple light scat- 
tering. 

The purpose of the present Letter is to fill up an im- 
portant gap in the quantum theory of multiple light scat- 
tering by developing a quantum description of dynamic 
multiple light scattering in fluctuating disordered media. 
In contrast to the previous work 0, S] , we take into ac- 
count the evolution of disordered medium (e.g., Brown- 
ian or more complex motion of scattering centers) in the 
course of a multiple scattering experiment. This leads 
to an interesting and reach problem of quantum dynam- 
ics, where the quantum fluctuations of light are coupled 
with the classical fluctuations of the medium. Our re- 
sults allow us to estimate the impact of disorder on the 
information capacity of quantum communication chan- 
nels and to lay a foundation for the use of nonclassical 
light as a probe of disordered media. 

We consider a slab of disordered medium that, for con- 
creteness, we assume to be a suspension of mobile, elasti- 
cally scattering particles (see Fig.[T]). For each frequency 
w, light incident on the slab can be decomposed in M 
transverse modes that we denote by indices a' and b' for 
modes on the left and on the right of the slab, respec- 
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FIG. 1: An ensemble of mobile scattering particles (circles) 
is illuminated from the left by a wave in the incoming mode 
a. All other incoming modes are in vacuum states (modes a' 
on the left and b' on the right). We are interested in the fluc- 
tuations of transmitted light in the outgoing mode b. Dotted 
line illustrates the multiple scattering with a mean free path 
I <^ L that light undergoes inside the medium. 



tively. All but one of the incoming modes are assumed 
to be in vacuum states. The state of the mode a' = a on 
the left can be arbitrary and it corresponds to the light 
that one sends into the medium in an experiment. We 
are interested in the properties of multiply scattered light 
in the transmitted, outgoing mode b. 

To quantize the electromagnetic fleld we use the same 
method as Beenakker et al. 7] and Lodahl et al. [§] 
who considered disordered media with immobile scatter- 
ing centers. An annihilation operator [uj) is associated 
with each mode and standard bosonic commutation re- 
lations are assumed 0, [1] : 



(1) 



where i and j run either over all incoming modes or over 
all outgoing modes. We now note that, as suggested by 
Fig. [Jl the fluctuating disordered medium acts as a beam 
splitter having M input and M output ports and time- 
dependent transmission and reflection coefficients Ta'b{t) 
and pb'b(t). By analogy with the work of Hussain et al. 
[3l who dealt with a 'regular' beam splitter (M = 2) 
having time-dependent properties, we assume that ra'&(t) 
and pb'bit) are frequency-independent and that they do 
not vary significantly during the typical time it takes for 
light to cross the disordered sample. The former assump- 
tion requires a sufficiently narrow bandwidth B of the in- 
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cident light, whereas the latter one is justified provided 
that the dynamics of the medium is slow. With these as- 
sumptions in hand, the input-output relations for aj(w) 
read 
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(2) 



where ai{t) denotes a Fourier transform of 0^(0;). The 
use of ai{t) is permissible provided that B is much 
smaller than the central frequency loq of the incident 
light [3, [l3 . Similarly to the case of immobile scatter- 
ers [71, y, Eq. ^ conserves commutation relations ((1]), 
which is due to the unitarity of the scattering matrix 
composed of Ta'bit) and pb'b{t) in the absence of absorp- 
tion. However, in contrast to the previous work 0, S], 
Eq. ((2]) implies a possibility of energy exchange between 
different spectral components upon transmission of light 
through a fluctuating disordered medium. This becomes 
clear when we back Fourier-transform Eq. ([2|) because 
products of time-dependent quantities become convolu- 
tions in the frequency space. Such an energy exchange 
is by no means a nonlinear phenomenon, but is a man- 
ifestation of Doppler effect due to scattering of light on 
moving particles. 

To establish a link between di(t) and measurable quan- 
tities, we define a flux operator ii{t) = al{t)ai{t) that, 
when integrated from t to t + T, yields a photon number 
operator hi{t^T) corresponding to the number of photo- 
counts measured by an ideal, fast photodetector illumi- 
nated by light in the mode i during some sampling time 
T [l3| • A relation between the average values of opera- 
tors hb{t,T) and ha{t,T) corresponding to the outgoing 
mode b and the incoming mode a, respectively, is read- 
ily found from Eqs. H]) and ([2]): {fib) = Tab(na)- Here 
Tab{t) = \Tab{t)\'^, the angular brackets denote the quan- 
tum mechanical expectation value, and the bar denotes 
averaging over disorder. By the ergodicity hypothesis, 
the simultaneous averaging (. . .) over both disorder and 
quantum fluctuations can be replaced by time averaging 
in a realistic experiment with a fluctuating disordered 
medium, provided that the incident light is stationary. 
On the contrary, there is no simple way of obtaining av- 
erages over either disorder or quantum fluctuations sep- 
arately. 

Although the relation between (hb) and {fia) is inde- 
pendent of the quantum state of the incident light, the 
latter is crucial for the second order statistics of the 
photon number operator. To illustrate this, we con- 
sider a particularly simple case of a stationary, quasi- 
monochromatic incident beam for which both {Ia{t)) and 
(: Ia(ti)Iait2) :), with : ... : denoting normal ordering 
[18|, can be considered time-independent on the scale of 
T. This includes a number of practically important ex- 
amples on which we will dwell in more detail: coherent 
light, continuous- wave quadrature squeezed light with 
a rectangular squeezing spectrum of width B <C 1/T 



around the central frequency of a perfectly monochro- 
matic mean field Aq exp{—iujot + i(p) [19|, light in a gener- 
alized continuum number state with a box-shaped 
envelope function $,{t) — (t,nax ~ imin)~^^^ exp(— iwoi), 
where imin ^ ^1,^2 < ^max ^21], and chaotic light with 
a long coherence time t^oh ^ T. Making use of Eqs. 
(HI) and ([2]), we find that the normalized variance of pho- 
ton number fluctuations in the outgoing mode b can be 
written as 
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Here Varfnf,] = {fi^) — {hb) , (ha) is the average photon 
number and Qa — Ya,T[na]/ (ha) — 1 is the Mandel pa- 
rameter [3l corresponding to the incident light beam. 
<5,2i_ = (2/T)Xf d^(l - t/T)Cf {t) is the normalized 
variance of time-integrated intensity in the mode b that 
one obtains from Eq. ^ by ignoring the quantum na- 
ture of light and dropping the operator 'hats' of a,i{t)] 

Cf{t) = Tab{ti)Tab{ti +t)/T^f,-l is the autocorrelation 
fimction of Tab{t). (^^i^gg represents the 'classical' part of 
6^ in the sense that it is determined uniquely by the (clas- 
sical) fluctuations of transmission coefficient Tab(t) due 
to the motion of scattering centers in the medium. 

The main conclusion following from Eq. ([3]) is that, in 
the general case, the variance of photon number fluctu- 
ations in transmission through a fluctuating disordered 
medium is not simply a sum of quantum and classical 
(due to disorder) contributions, but rather a much more 
complicated object. In particular, the last term of Eq. 

mixes classical fluctuations and quantum noise in a 
multiplicative way. Only for incident light in the coher- 
ent state Qa = and Eq. ([3]) reduces to a sum of shot 
noise l/{hb) and classical noise S'^^^^^. This is the mini- 
mum value of 6'^ that can be obtained for light in a state 
admitting classical description. An 'excess noise' ^excess 

— defined as a difference between Sf and l/ihb) + <5^iass 

— arises for other states. For chaotic light, for example, 
Qa = {na) and (J^xcess = 1 + (^ciass- ^ote that for both 
coherent and chaotic states we do not really need the 
quantum model developed here, but can instead use the 
semiclassical Mandel's formula as it has been done 
in Ref. [12]. 

We now turn to nonclassical light — light in states 
that cannot be described classically and to which the 
Mandel's formula does not apply. For incident light in 
the number state \na,£,), Qa = ~1 and the excess noise 
is negative: the fluctuations of transmitted photon num- 
ber are suppressed below the simple sum of shot noise 
l/{hb) and classical fluctuation S'^i^^^. 5l becomes even 
less than the bare shot noise l/ihb) if < 1 + ^ciLs- 
interesting situation occurs for the single-photon number 



FIG. 2: Excess noise of photon number hb in the transmitted 
mode b for the incident light (mode a) in the continuous- wave 
amplitude-squeezed state with a narrow rectangular squeez- 
ing spectrum of width B <g 1/T. The excess noise is defined 
as a difference between the variance of photon number fluctu- 
ations the shot noise l/lfib) and the classical noise (Jciass- 
It is shown as a function of sampling time T and squeez- 
ing parameter r for a suspension of small dielectric particles 
in Brownian motion (diffusion coefficient Db) ss depicted in 
Fig. [T] The unit of time is to = /GPDbL^ , where / is the 
mean free path, L ^ Hs the sample thickness, and k ^ l/l is 
the wave number of light. The intensity carried by the mean 
field Aq is set equal to B. Thick black lines show the excess 
noise as a function of sampling time for r — 0, 0.2, 0.4, 0.6, 
0.8 and I.O. 



state = 1): the relative fluctuation of transmitted 
photon number becomes independent of iJ^iass- Another 
important example of nonclassical light is the continuous- 
wave squeezed Hght for which Qa = {ha)f{r, (j), A'^/B, 1) 
[l^ . where r is the squeezing parameter and we define 
an auxiliary function f{r,(j>,x,y) — [xy{e~'^'^ cos^ (j) + 
e^^ sia^ (f) — 1) -I- sinh^ r cosh 27-]/(a; + sinh^r)^. For 
the amplitude-squeezed light (0 = 0) the Mandel pa- 
rameter is negative if r < tq, with rg depending on 
A^/B, leading to (Joxccss < 0' similarly to the num- 
ber state. This is illustrated in Fig. [2] for the case of 



C^{t) = (t/to) / sinh yJtJtQ which corresponds to a sus- 



pension of small Brownian particles [Bl, lij (the character- 
istic decay time to is defined in the caption of Fig. [2]). 
With A^/B — 1 chosen for this figure the excess noise is 
negative for r < ro ~ 0.6. 

Equation ^ can be used to estimate the effect of dis- 
order on the information capacity of quantum commu- 
nication channels operating in disordered environments. 
As an example, let us consider a communication channel 
connecting Alice on the left of the disordered medium in 
Fig. [T] and Bob on the right. Alice communicates with 
Bob by sending a sequence of number states, each of du- 



FIG. 3: Autocorrelation function of photon number fluctu- 
ations in transmission of amplitude-squeezed light through a 
suspension of small dielectric particles in Brownian motion for 
the same values of parameters as in Fig. [2] and to = 0.1/ B. 
The dashed red line illustrates that C^(t) is independent of r 
for t = l/B. 



ration T and containing up to n^^^ photons, in the in- 
coming mode a. Bob reconstructs the message of Alice by 
counting photons in the outgoing mode h during the same 
time intervals T. By using the theory of linear bosonic 
communication channels [20'| and Eq. ([3]), we find that 
the information capacity of such a channel is 
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where we drop higher-order terms in n^^^^Tab ^ 1, which 
is justified by the condition Tab 1 typical for transmis- 
sion of light through a diffusely scattering medium. As 
could be expected, the (classical) fluctuations S^^^^ re- 
duce the capacity of the channel. However, this reduction 
is only second-order in n^^^Tab ^ 1- The drop of capac- 
ity in the presence of disorder is therefore mainly due to 
the reduced average value of transmission coefficient Tab 
rather than to the fluctuations of the latter. Interest- 
ingly, C appears to be independent of 5^^^^^ if n^'^^ = 1 
(Alice sends either f or no photons). 

Correlation of photon numbers in the same mode b but 
in two different, nonoverlapping time intervals, 



{hi,{ti,T)nb{ti+t,T)) 
{hb{ti,T)) {hb{h+t,T)) 



1, 
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is also affected by the quantum state of incident light and 
can be reduced to 



(6) 
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where C°(i) is the autocorrelation function of photon 
number fluctuations in the incident mode a. In order 
to put Eq. © in a simple form, the typical decay time 
of C^{t) is assumed to be much longer than T. Equa- 
tion ([5]) is more general than Eq. ([3]) and requires neither 
stationarity nor monochromaticity of incident light. For 
light in the coherent state C^{t) = and Cl{t) = Cf{t). 
This result actually lies at the heart of diffusing-wave 
spectroscopy (DWS) and other photon correlation tech- 
niques that identify the measured C^{t) with C^{t) and 
interpret it in terms of classical wave scattering. 

For nonclassical states of incident light C^{t) is not 
equal to C^{t). For the number state |na,Oj ^ni^) = 
— l/ua and C^{t) can take negative values and becomes 
independent of time for ria — 1. These strange proper- 
ties are direct consequences of the conservation of total 
number of photons — a conservation law which does not 
apply to the coherent state where the number of pho- 
tons is not a good quantum number. For squeezed light 
we find C^{t) = f[r,(j),Al/B,sinc{-KBt)]. As illustrated 
in Fig. [31 squeezing induces oscillations of C^^{t) and de- 
creases its value at short t. The latter reaches a minimum 
at some r = ri (ri ~ 0.3 in Fig. [3]) and then grows again, 
exceeding 1 at large r. The results presented in Fig. [3] 
open a way to performing DWS with squeezed light; the 
necessary experimental setup is completely analogous to 
the 'traditional' DWS one [1, @|, except that a source of 
squeezed light should be used instead of a laser. 

Provided that Tab{t) is replaced by Ta[t) = J2bTab(t) 
in the definition of C^''{t), the main results of this Let- 
ter — Eqs. ([3]) and ^ — also hold for the variance and 
the autocorrelation function of fi = fib — a photon 
number operator corresponding to the 'total transmis- 
sion' measurement. Equations ([3]) and ([6]) thus implic- 
itly account for all possible long-range (in space or in 
time) correlations of Tab{t) that were extensively stud- 
ied previously [^, [l3| ■ The effect of long-range correla- 
tions on photon number fluctuations in total transmission 
through a suspension of scattering particles in Brownian 
motion has been recently studied by Balog et al. 12 1 



for incident light in the coherent state, and the result 
following for this case from Eq. ([3]) has been confirmed 
experimentally. 

In conclusion, we developed a quantum description of 
dynamic multiple light scattering in fluctuating disor- 
dered media. We calculated the fluctuation and the auto- 
correlation function of photon number operator in trans- 
mission through a slab of mobile scattering particles. We 
applied our results to estimate the effect of disorder on 
the information capacity of a quantum communication 
channel operating in a disordered environment and to dis- 
cuss the possibilities of performing diffusing-wave spec- 



troscopy with squeezed light. Our results suggest that 
using nonclassical light (single photons, light in number 
states, squeezed light) opens interesting perspectives in 
optics of disordered media. 
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